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Abstract 

NOON states, states between two modes of light of the form |A^, 0) + 6**^10, A^) allow for super- 
resolution interformetry. We show how NOON states can be efficiently produced in circuit quntum 
electrodynamics using superconducting phase qubits and resonators. We propose a protocol where 
only one interaction between the two modes is required, creating all the necessary entanglement 
at the start of the procedure. This protocol makes active use of the first three states of the phase 
qubits. Additionally, we show how to efficiently verify the success of such an experiment, even 
for large NOON states, using randomly sampled measurements and semidefinite programming 
techniques. 
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INTRODUCTION 



The particle-wave duality of quantum physics has striking consequences on the physics of 
electromagnetic radiation, necessitating a photonics description of electromagnetic radiation, 
spanning frequencies from 7 quanta via visible light to microwaves fields. The phenomena 
of classical optics are fully characterized by only a tiny fraction of the wealth of states 
attainable in the quantum world. Creating and quantifying these others, the nonclassical 
states of light, has defined the field of quantum optics [1, 2]. 

Recently, this quantum-optical paradigm has been extended to the field of circuit quantum 
electrodynamics (cQED) [•>, 1] where nonlinear superconducting circuits [5] are used to 
manipulate the quantum state of microwaves in waveguides in much the same way as atoms 
are used to manipulate optical fields in traditional quantum optics. Unlike atoms, however, 
superconducting circuits can be tuned over a broader range of parameters, allowing for 
highly flexible control. The current state of the art allows for the full control of quantum 
states up to 12 photons [(>] as well as violations of locality in a Bell experiment [7, 8]. Small 
quantum information algorithms have also been successfully demonstrated [9]. 

A new frontier in this program is the realization of nonclassical states that are entangled 
between modes in distinct waveguides. A prominent example is the NOON state [10-12] 

|^) = i=(|iV,0) + e^^|0,iV)), (1) 

i.e., an equal coherent superposition of photons in the first mode and vacuum in the 
second with the reverse, vacuum in the first and A^ photons in the second. NOON states 
were introduced in quantum metrology as a means of surpassing the shot noise limit for 
phase sensitivity achievable by classical light by a factor of 1/ \/N, which instead saturates 
the Heisenberg bound. At optical frequencies, NOON states up to A^ = 4 have been realized 
to date, relying on cascaded pair sources and postselection [13, 14]. Schrodinger cat states 
interpreted as NOON states up to A^ = 10 have been realized in spin systems, along with 
an enhanced sensitivity to the magnetic field [15]. 

We propose a scheme to create NOON sates in cQED. In cQED, sources that create pairs 
of photons independent of the number of photons already in the mode are currently not 
available. These would be needed for direct implementation of schemes developed at optical 
frequencies. We will circumvent this restriction by utilizing an additional energy level in 
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the superconducting qubit we are considering for our artificial atom. We also place bounds 
on the size of NOON state that can be achieved and propose a scheme how to verify the 
preparation of a NOON state using single-qubit measurements. 

SYSTEM AND HAMILTONIAN 

The system under consideration consists of two microwave cavities, each coupled to its 
respective phase qubit. The qubits are also coupled to each other, either directly [16-18] or 
by a coupling cavity [19, 20]. This being a non-essential difference, we are going to focus 
on the case of cavity-mediated interaction, where by programming the frequency in phase 
qubits one achieves tuneability . A schematic this system is given in Fig. 1, and is described 
by the Hamiltonian 

H{t) = Ho + H^Q.it) + H^Q,it) + H^ciit) + H^t)- (2) 
The static (or 'drift') part of the Hamiltonian can be decomposed as 

Ho=^ Hqi + ^ Hcj + Hqicj. (3) 

i=l,2 i=l,2,C i,j 

Here, Hqi = En |1). (l|. + £'2i |2)j (2|. is the Hamiltonian for the phase qubit truncated to its 
three lowest energy levels |0), |1), and |2). Higher states will never be resonant during the 
whole protocol. The relevant cavity modes are described by Hcj = ^j<^]<^j with j = 1,2, C 
and Oj and aj being the photon raising and lowering operators respectively. With the 
same index space, the qubit-cavity couplings can be written in rotating wave approximation 
(RWA), assuming g^^^ <^ uj. En as 

HQ^c, = {gff (0|, + |2), (1|,) a, + h.c. (4) 

which is a straightforward three-level generalization of the Jaynes-Cummings (JC) interac- 
tion [21]. 

The system is primarily manipulated by controlling the phase qubits. Driving 
one of the qubits with an external field results in a control of the form, Hq^it) = 
(|0)j (l|j + A |l)j (2|^ -|- h.c). Here, A ~ \/2 is a constant whose precise value depends on 
details of the qubit. Generally, Q{t) will be a resonant microwave drive Q{t) = Qq cos{ujdt) . 
Additionally, changing the bias current gives a time-dependent En and E2i approximately 
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FIG. 1: Schematic of an experimental setup for generating NOON states consisting of two phase 
qubits, a coupling cavity, and two data cavities. The NOON state will be between the data cavities, 
which in the end are disentangled from the rest of the system. The boxed x's label Josephson 
junctions and the small circles external ports to which driving fields can be applied 

following Eli = pE2i with a fixed p < 2. Changing the energy of the qubit will primar- 
ily be used as a switch to move interactions in and out of resonance. Finally, the cavities 
themselves can be driven externally by microwave sources H^- = Ai(t)(a| + di). 



NOON STATE PREPARATION 
Coupling and decoupling 

In order to understand the coupling Hamiltonian, it is helpful to transform to the inter- 
action representation with respect to the Hqi and Hcj- In that frame 

^.c, = {gff (0|.e<* + |2), (l|,e<') a, +h.c.. (5) 

with the detuning aJ^^^ = El''\t)-El''~^\t)-Uj. If A^f (t) > g-j^ y/rij + 1, a standard rotat- 
ing wave argument allows us to neglect the terms that couple \k)^ \n)j ^ \k — 1). \n + 
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The implication is that by rapidly tuning , one can switch the coupling between the 
cavity and the qubit on or off. This has been successfully demonstrated in a multitude of 
experiments [6, 19, 22]. The two regimes of interest for the following protocol are the case 
when both detunings are large, leading to uncoupled systems, and the case when A[f = 
and Aj-^'' ^ gf^ \f2, which reduces to the standard resonant JC interaction between states 
|0). , and the qubit only, see fig. 2, lower right. 

^ic,=^?ll^(|l).(0|.a, + |0),(lU. (6) 



Bell pairs 

Preparation of Bell states in this setting has already been demonstrated in Refs. [7, 8]. 
It works as outlined in fig. 2, lower left using the notation l^i, ?^c, foi' qubit 1, cavity 
C, and qubit 2. Starting from all elements decoupled by detuning, we let the system relax 
into its ground state |0,0,0). A resonant vr-pulse on generates the state |1,0,0). Then, 
we bring qubit 1 into resonance with the coupling cavity, A^)^ = 0, for a time t\ such that 
g^y^Tijl = IT /2 to create (|1,0,0) — i |0, l,0))/-\/2. Next, we move qubit 1 out of resonance 
and qubit 2 into resonance for ^2C^2/2 = vr leading to the state (|1,0,0) + z |0, 0, l))/y2, 
which is a Bell state between qubits 1 and 2, decoupled and disentangled from the cavity, 
see fig. 2. After this step, there is no requirement to couple the qubits again throughout 
the whole protocol, and so we have no further use for the coupling cavity. All that remains 
is to separately generate a map on each subsystem, (qubit+cavity)^ denoted by \qi,ni), that 
transforms |0,0)^ — ?■ |0,A^)^ and |1,0)^ — |0,0)^. All necessary entanglement has already 
been created, since an arbitrary NOON state has exactly the same amount of entanglement 
as a Bell pair. This is entanglement between the subsystems (i.e. the two cavity/qubit 
subsystems), as opposed to an entanglement measure that treats the photons themselves 
as particles that could be made distinguishable by later splitting them into different modes 
[10, 11]. 

Shelving 

In order to selectively increase the photon number of the qubit ground-state branch of 
the superposition while keeping the excited-qubit branch in vacuum, we shelve the excited 
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FIG. 2: The NOON-state preparation protocol. Top: Overview over the steps, highhghting qubits 
1 and 2 and cavitites 1, 2 and C. Parties 1 and 2 only interact during the initial Bell state 
preparation. Bottom left: Bell-state preparation by interacting qubits 1 and 2 with the cavity. 
Qubit 1 has already undergone a vr-rotation and is in the state |0) + |1) /y/2. The time axis shows 



interation times ri and T2 necessary for performing a ViSWAP operation with the first and an 
iSWAP with the second qubit. The frequency axis indicates that the interaction is controlled by 
moving the qubits in and out of resonance with the coupling cavity. Bottom right: Energy scheme 
indicating the cavity energy ladders for the three different phase qubit states, assuming resonance 
between the cavity and the -H- 1 transition of the qubit. In the last 'unshelving' step of the 
protocol, |0,n) and — 1) are swapped. 
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state 1 1,0)^ in the second excited state |2,0). by selectively driving the qubits at frequency 
hwd = E2^i — Ei^i. The |2,0),- state is considered 'shelved' in the sense that the Jaynes- 
Cummings (JC) interaction that is resonant with the transition |0, 1) ■ ^ |1, 0)^ is far detuned 
form the transition |1, 1)^ ^ |2, O),- as long as £'2,* — 2_Ei j ^ gl^', see fig. 2. The |2, 0)^ state 
is also left invariant by rotations on the first two qubit levels assuming that the nonlinearity 
leaves it far off-resonant, or through the use of pulse shaping techniques such as DRAG [23]. 

Climbing the Jaynes-Cummings ladder 

As a next step, the non-shelved state |0, 0)^ is transformed into the state |0, — 1)^. This 
is a special case of the Eberly-Law protocol [21] that has been demonstrated in refs [22, 25]. 
Each step from \Q,n — 1)^ to |0,n)^ is realized by first exciting the qubit by a resonant vr- 
pulse, generating |1, n — 1).. Next, the qubit is moved into resonance with its cavity for time 
gf^ ^fntjl = IT, moving the excitation from the qubit to the resonator, see fig. 2. This is 
repeated until |0, A^ — 1)^ is reached. At the same time, due to efficient shelving, the other 
branch of the superposition stays in |2,0)^ as the |2,0)- ^ transition is off-resonant 

with the cavity, see fig. 2, lower right. 

Unshelving 

The procedure so far brings the total state, written in notation |ni, 712, qi, q^-, nc) in state 
(|A^ - 1, 0, 0, 2, 0) - |0, A^ - 1, 2, 0, 0))/y2. At this point it may be counter-intuitive as to 
why we stopped at A^ — 1 photons in the previous section, but the simple explantation is 
that to unentangle the qubits from the cavities we must use the JC-interaction which will 
inject an additional photon into our state. The unshelving procedure is as follows. Applying 
a vr-pulse to both qubits on the \2)- ^ transition, followed by vr-pulses on the transitions 
o |0)j permutes the qubit states such that \2)- |0)j and |0). — )■ This results in 
the state (|Ar - 1, 0, 1, 0, 0) - |0, A^ - 1, 0, 1, 0))/V2. 
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Disentangling 



Now we couple both qubits' to their respective cavities, A^^^ = 0, for time gll''y/Nt/2 



= 71. 



If both the cavity and the qubit are in their ground states they are left invariant by the JC 
interaction. On the other hand, the branch of the wavefunction that contains |1,A^ — 1)^ 
becomes |0, A^)^. This step disentangles the cavities from the qubit and leads to the desired 



IMPACT OF DECOHERENCE, MAXIMAL SIZE ATTAINABLE 

While there is no principal bound on the maximum size N of NOON state this protocol 
can achieve, there are practical limitations. One will be related to the precision of the 
control pulses that can be applied, which in principle can be very high using pulse shaping 
[23, 26, 27]. More fundamentally, decoherence limits the size of the state attainable. There 
are two crucial limitations: One is that the qubit state has to be shelved in the second 
excited state for a long time and the other is that the superposition should not dephase. 
Given that his happens if any of the two qubits decays, we can go up to T2/2. The total 
duration of our protocol can be simply calculated as. 



where T^abi is the time for a single qubit Rabi oscillation. Assuming T2 ^ 200ns [28] and 
Tj^abi — 10ns and 1/g — 20ns, this allows NOON states up to = 4. These numbers fall 
well below the number of photons already stored and investigated in such resonators, so 
photon loss will not be a limiting factor [6, 29]. In transmons, which also offer an accessible 
third level [30], coherence times are longer [31] and larger NOON states may be attainable. 

STATE VERIFICATION 

Decoherence and imperfect control ensure that in the laboratory setting, the prepared 
state described by a density matrix pcxp is not a perfect NOON state. The fidelity between 



NOON-state 



|z^) = 4^(|iV,0)-|0,iV))®|0,0,0). 



(7) 




(8) 
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the desired and actual state can be obtained by measuring the overlap 
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(9) 



This is not a trivial task with the measurement tools at our disposal, which consist of driving 
the system through some simple dynamics followed by qubit measurements. In this section, 
we will show that we can create a set of measurement operators that is informationally 
complete, that is they span the operator space 5\i{d) where d is the dimension of the problem 
under consideration. Here d = {N + 1)^ when we consider just the microwave cavities. The 
outcomes from an informationally complete set of measurements can be used to perform 
tomography on the quantum state, but we will also present a method of bounding the 
fidelity directly. 

Restricting our attention to a single cavity, the types of measurements we consider consist 
of driving the cavity by a coherent displacement a, switching on the JC Hamiltonian for a 
time gf^tjl = r, and then measuring the qubit in the 2;-basis by measuring the switching 
probability out of state |1) [5, 28]. When we make the simplifying assumption that the qubit 
is initially uncoupled from the cavity and in its ground state, we can express this observable 
on the microwave cavity as 



By sampling different values of the dimensionless time r, with the amplitude set to a = 0, 
one can only infer the diagonal elements of the density operator {n\ p \n). The injection of 
a is crucial because diagonal measurements alone do not distinguish a NOON state from 
the mixed state p = |(|A^, 0) (A^, 0| + |0,A^) (0,A^|), the latter being useless for quantum 
metrology. In [6], these types of measurements were used to construct the Wigner function 
of the cavity, which is a complete description of the quantum state of the cavity mode. 
This implies, that at least to some finite truncation of the cavity, the set of measurements 
Ai{a,T) are informationally complete in that one can can construct a set of (A^ + 1)^ — 1 
linearly independent observables spanning the space 5u{N +1). To generate such a set it is 
sufficient to simply sample pairs {aj,Tj) from some random distribution. This will almost 
surely generate a linearly independent set of measurement operators. It is not impossible 
that there are optimized choices of operators that are close to orthonormal, thus optimizing 
the sampling of phase space, but in practice it's most often easier to instead just oversample 
from our distribution on (a, r). 




(10) 



n 
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Correlating the results of measurements from the two individual cavities yields the set 
of operators 7^2(0, t, a', r') = J\4{a, r) (g) Ai{a', r'). This is a complete set of operators on 
the joint space, however, since the Hilbert space of the two cavity system has dimension 
(A^ + 1)^, reconstructing the entire density matrix requires measuring at least (A^ + 1)^-1 
of these operators. To perform full tomography we repeat each measurement a number of 
times to obtain measurement statistics 

Mj = Tr (>l2(a„ r,, a'j, rppex-p) + cXjWj, (11) 

where the final term is a Gaussian white noise variable with variance aj. In principle, one 
can arbitrarily decrease the strength of this noise by repeated measurements, however, it 
is generally more efficient to fix aj and instead sample more measurement operators from 
our distribution on Ai{aj,Tj,a'j,Tj). From such a measurement record one can obtain a 
maximum likelihood estimate of the quantum state through a least squares fit [32]. The 
resulting state may have unphysical, negative eigenvalues which need to be compensated 
for by either simply setting them to zero or by using slightly more involved semi-definite 
programming techniques [-33, 34] to find the physical state that is most consistent with 
the measurement outcomes. The assumption that the qubits were initially in their ground 
states was not essential and can be circumvented by adding to our set of measurements by 
additionally performing some unitary operations to the superconducting qubit immediately 
prior to measurement. This increases the number of independent measurements necessary 
for tomography to 16(A^ + — 1 for qubits, or 81(A^ + 1)^ — 1 when we also account for 
the second qubit excited state. 

Placing bounds on the experimental value of J-" can be accomplished with a smaller set 
of measurement outcomes. We can lower bound J-", defined in eq. (9), by directly solving 
the optimization problem 

minj^ (12) 

p 

subject to 

Mj = Tr (^M{aj, t„ a'^, t'.)p^ , p > 0. (13) 

where we have neglected the noise on the measurements. The second constraint ensures 
that p is a physical density matrix, and thus we will find the state with the lowest overlap 
with the NOON state that nevertheless agrees with the outcomes from our measurements. 
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This optimization is of the form of a semidefinite program and can be solved very efficiently 
numerically [33]. We used the matlab program yalmip running the solver SeDuMi. Addition- 
ally, since the objective of this optimization (J-", eq. 9) is linear, solving the maximization 
problem (i.e. finding an upper bound), is of the same difficulty. 

To determine the utility of this bounding scheme we ran some simple numerical sim- 
ulations, the results of which are displayed in Fig. 3. We place upper and lower bounds 
on states that are convex sums of the NOON state and the maximally mixed state, or 
Pcxp = PPnoon + (1 — For a given overlap with the NOON state, these mixtures are 

generally the states that are hardest to bound, and so represent a worst case test of the 
protocol. The measurements were chosen by sampling the Tj uniformly from the interval 
[0, tt] and the a's uniformly from the interior of the circle on phase space with radius A^. 
For a simple test we looked at NOON states with photon numbers going up to = 6, 
where even the longest semidefinite optimizations were solved in about a minute on a laptop 
running matlab. 

In the case of randomly sampled measurements, the lower bound has a rapid initial rise 
with respect to the fraction of the possible independent operators measured. The time 
for convergence of the upper and lower bounds is highly dependent on the fidelity of the 
experimental state, with high fidelity states converging very quickly. In order to get an fairly 
exact estimate of the fidelity for a density matrix far from the desired state, this method 
may not be much more efficient than full tomography, however, it is often useful to quickly 
generate bounds out of small data sets in order to diagnose state generation protocols. Also, 
we suspect that this protocol is more amenable to optimizing the measurement operators 
than full tomography, and so it should be possible to get faster convergence than can be 
achieved by random sampling. 

CONCLUSIONS 

In conclusion, we have shown how to determinisitally create large NOON states of mi- 
crowave photons in superconducting resonators, and how to verify the success of this ex- 
periment in a more efficient way than full state tomography. It can be shown, that this 
reconstruction protocol can be extended to account for detector noise. 

NOON states in the optical domain have been proposed for super-resolution in quantum 
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FIG. 3: Numerical analysis of bounding the fidelity of NOON state construction with an incomplete 
set of measurement outcomes. All simulations are for noiseless measurements and assume the qubits 
are initially in their ground states, (a) Upper and lower bounds for the overlap between the 4- 
photon NOON state and various mixtures of the this NOON state and the maximally mixed state, 
plotted vs. the number of randomly sampled measurements performed, (b) Lower bounds on the 
NOON state fidelity for mixtures of different dimension. The x-axis of this plot is the number of 
measurement results divided by the dimension of 5u{d) for the different NOON states, (c) The 
difference between the upper and lower bounds for mixtures of the 6-photon NOON state and the 
maximally mixed state plotted versus the number of measured observables. (d) The number of 
measurements necessary to give a separation of 5% between the upper and lower bounds, plotted 
versus the actual fidelity of the experimental state. 
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lithography and interferometry [10, 11, 15]. With wavelengths of microwave photons in our 
case being on the order of centimetres, lithograqhy is unlikely for microwave photons. Mi- 
crowave interferometry has applications in plasma diagnostics [35] and astrophysics relating 
to the cosmic microwave background. The application of NOON states to these systems 
remains speculative, and will be explored in future work. 

We acknowledge useful discussions with J.M. Martinis and A. Cleland, who indepen- 
dently developed a similar protocol for NOON state creation, as well as discussions with M. 
Mariantoni, H. Wang, F. Motzoi, and B. McNamara. Work supported by NSERC through 
Quantum Works and the discovery grants program and, in parts, by DARPA through the 
QuEST program and lARPA through the MQCO program. 



[1] R. Glauber, Phys. Rev. 131, 2766 (1963). 

[2] K. Cahill and R. Glauber, Phys. Rev. 177, 1882 (1969). 

[3] A. Blais, R.-S. Huang, A. Wallraff, S. Girvin, and R. Schoelkopf, Phys. Rev. A 69, 062320 
(2004). 

[4] R. Schoelkopf and S. Girvin, Nature 451, 664 (2008). 
[5] J. Clarke and F. Wilhelm, Nature 453, 1031 (2008). 

[6] M. Hofheinz, H. Wang, M. Ansmann, R. C. Bialczak, E. Lucero, M. Neeley, A. D. O'Connell, 

D. Sank, J. Wenner, J. M. Martinis, et al.. Nature 459, 546 (2009). 
[7] M. Ansmann, H. Wang, R. C. Bialczak, M. Hofheinz, E. Lucero, M. Neeley, A. D. O'Connell, 

D. Sank, M. Weides, J. Wenner, et al, Nature 461, 504 (2009). 
[8] J. Chow, L. DiCarlo, J. Gambetta, A. Nunnenkamp, L. Bishop, L. Frunzio, M. Devoret, 

S. Girvin, and R. Schoelkopf, arXiv:0908.1955. 
[9] L. DiCarlo, J. Chow, J. Gambetta, L. Bishop, B. Johnson, D. Schuster, J. Majer, A. Blais, 

L. Frunzio, S. Girvin, et al, Nature 460, 240 (2009). 
[10] A. Boto, P. Kok, D. Abrams, S. Brownstein, C. Williams, and J. Dowling, Phys. Rev. Lett. 

85, 2733 (2000). 

[11] P. Kok, H. Lee, and J. Dowling, Phys. Rev. A 65, 052104 (2002). 
[12] B. Sanders, Phys. Rev. A 40, 2417 (1989). 

[13] M. Mitchell, J. Lundeen, and A. Steinberg, Nature 429, 161 (2004). 



13 



[14] T. Nagata, R. Okamoto, J. O'Brien, K. Sasaki, and S. Takeuchi, Science 316, 726 (2006). 
[15] J. Jones, S. Karlen, J. Fitzsimons, A. Adravan, S. Benjamin, G. Briggs, and J. Morton, Science 
324, 1166 (2009). 

[16] A. Berkley, H. Xu, R. Ramos, M. Gubrud, F. Strauch, P. Johnson, J. Anderson, A. Dragt, 

C. Lobb, and F. Wellstood, Science 300, 1548 (2003). 

[17] R. McDermott, R. Simmonds, M. StefFen, K. Cooper, K. Cicak, K. Osborn, S. Oh, D. Pappas, 

and J. Martinis, Science 307, 1299 (2005). 
[18] M. StefFen, M. Ansmann, R. C. Bialczak, N. Katz, E. Lucero, R. McDermott, M. Neeley, 

E. M. Weig, A. N. Cleland, and J. Martinis, Science 313, 1423 (2006). 
[19] M. Sillanpaa, J. Park, and R. Simmonds, Nature 449, 438 (2007). 

[20] J. Majer, J. Chow, J. Gambetta, J. Koch, B. Johnson, J. Schreier, L. Frunzio, D. Schuster, 

A. Houck, A. Wallraff, et al.. Nature 449, 443 (2007). 
[21] E. Jaynes and F. Cummings, Proc. IEEE 51, 89 (1963). 

[22] M. Hofheinz, E. Weig, M. Ansmann, R. Bialczak, E. Lucero, M. Neeley, A. O'Connell, 

H. Wang, and J. Martinis, Nature 454, 310 (2008). 
[23] F. Motzoi, J. Gambetta, P. Rebentrost, and F. Wilhelm, Phys. Rev. Lett 103, 110501 (2009). 
[24] C. K. Law and J. H. Eberly, Phys. Rev. Lett. 76, 1055 (1996). 

[25] J. Fink, M. Goppl, M. Baur, R. Bianchetti, P. Leek, A. Blais, and A. Walhaff, Nature 454, 
315 (2008). 

[26] N. Khaneja, T. Reiss, C. Kehlet, T. Schulte-Herbriiggen, and S. Glaser, J. Magn. Reson. 172, 
296 (2005). 

[27] A. Sporl, T. Schulte-Herbrueggen, S. Glaser, V. Bergholm, M. Storcz, J. Ferber, and F. Wil- 
helm, Phys. Rev. A 75, 012302 (2007). 
[28] J. Martinis, Quantum Information Processing 8, 81 (2009). 

[29] H. Wang, M. Hofheinz, M. Ansmann, R. Bialczak, E. Lucero, M. Neeley, A. O'Connell, 

D. Sank, M. Weides, J. Wenner, et al., Phys. Rev. Lett. 103 (2009). 

[30] J. Koch, T. Yu, J. Gambetta, A. Houck, D. Schuster, J. Majer, A. Blais, M. Devoret, S. Girvin, 

and R. Schoelkopf, Phys. Rev. A 76, 042319 (2007). 
[31] A. Houck, J. Koch, M. Devoret, S. Girvin, and R. Schoelkopf, Quant. Inf. Proc. 8, 105 (2009). 
[32] M. Paris and J. Rehacek, eds.. Quantum State Estimation (Springer, 2004). 
[33] L. Vandenberghe and S. Boyd, SIAM Review 38, 49 (1996). 

14 



[34] S. T. Merkel, C. A. Riofrio, S. T. Flammia, and I. H. Deutsch, Phys. Rev. A 81, 032126 
(2010). 

[35] M. Laroussi, Int. J. Infrared and Millimeter Waves 8, 1501 (1999). 



15 



